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on  two  important  aspects  of  shaped  charge  mechanics,  jet  breakup  and  broken  jet 
penetration.  The  jet  breakup  problem  is  studied  from  a stability  approach. 

A one-dimensional  stability  theory  for  stretching  plastic  jets  is  presented. 
This  theory  yields  a stability  criterion  which  depends  on  the  stress-strain 
curve  of  the  material.  It  also  gives  a critical  disturbance  wavelength; 
disturbances  of  this  wavelength  grow  faster  than  all  others.  In  addition,  the 
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^Stretching  rate  and  increases  with  increasing  flow  stress  to  density  ratio, 
Comparison  with  two-dimensional  finite  difference  calculations  indicates  the 
general  theoretical  trends  are  correct. 

The  second  study  examines  the  problem  of  broken  jet  penetration  by 
using  two-dimensional  calculations  to  observe  the  effects  of  sequential 
projectile  impact  on  the  bottom  of  a hole  in  a target.  Past  experiments  show 
decrease  in  penetration  caused  by  the  increase  in  gap  between  projectiles. 
Calculations  of  lead  projectiles  upon  lead  targets  do  not  show  the  interference 
phenomenon  which  is  expected  to  account  for  the  decrease  in  penetration. 
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I. 


INTRODUCTION 


For  the  past  several  years,  considerable  effort  has  been  placed  on 
determining  the  basic  mechanisms  in  shaped  charge  Jet  formation,  breakup, 
and  penetration.  Various  approaches,  including  experimental,  analytical, 
and  numerical,  have  been  used  to  treat  the  shaped  charge  problem.  Several 
aspects  of  shaped  charge  mechanics,  especially  the  breakup  phenomenon, 
are  still  not  well  understood.  The  present  study  will  focus  on  two  of 
these  aspects,  jet  breakup  and  broken  jet  penetration. 

In  [1-3],  we  have  presented  the  results  of  initial  studies  on  the 
jet  breakup  phenomenon.  Formulas  for  jet  strain  and  jet  radius  were  de- 
veloped and  radiographs  of  broken  jets  measured.  This  work  was  extended 
in  [4]  to  include  a larger  cross-section  of  experimental  data  which  yielded 
an  empirical  correlation  (design  curve)  indicating  that  jet  breakup  time 
increases  with  increasing  jet  radius,  a general  trend  expected  from  linear 
scaling.  Further,  because  of  the  resemblance  of  shaped  charge  jet  breakup 
and  liquid  jet  breakup,  the  breakup  problem  was  also  approached  as  an  in- 
stability problem.  In  [4],  we  have  shown  that  the  classical  stability 
analyses  do  not  apply  to  shaped  charge  jets.  We  then  performed  two-dimen- 
sional finite-difference  calculations  using  the  HEMP  code  [5]  to  simulate 
stretching  elastlc-perf ectly  plastic  jets.  These  calculations  [4]  show 
that  the  shaped  charge  jet  is  indeed  subject  to  an  instability  which  causes 
surface  disturbance  growth  and  eventual  Jet  breakup.  Furthermore,  it  was 
found  that  the  driving  force  for  the  instability  was  caused  by  the  material 
strength  of  the  jet  and  not  the  surface  tension  as  in  a classical  liquid 
jet.  The  existence  of  a critical  (i.e.  most  unstable)  wavelength  in  shaped 
charge  jets  was  indicated  by  this  numerical  study.  Related  work  by  Karpp 
and  Simon  [6]  has  also  shown  the  existence  of  a critical  wavelength  and  the 
importance  of  jet  strength.  These  studies,  [4]  and  [6],  have  shown  the 
range  of  the  critical  wavelength  to  jet  diameter  ratio  at  its  time  of 

initiation  to  be  1.5  < (X/d)  . <2.5.  Of  course,  by  the  time  the  jet 

crit 

breaks,  this  surface  "wave"  will  be  stretched  out  so  that  the  final  aspect 
ratio  of  the  particles  will  be  larger. 

In  this  report  we  will  present  the  results  of  our  current  study  on 
jet  breakup.  This  study  is  a continuation  of  the  stability  analysis  of 
shaped  charge  jets.  A one-dimensional  stability  theory  for  stretching 
plastic  jets  is  presented  which  includes  the  essential  effects  controlling 
stability  indicated  by  the  two-dimensional  numerical  computations.  Further, 
two-dimensional  calculations  of  additional  factors  affecting  jet  stability 
are  presented. 

The  second  topic  of  this  report,  sequential  projectile  impact  and 
penetration,  is  an  idealized  model  which  simulates  broken  shaped  charge  jet 
penetration.  In  the  past,  experiments  were  conducted  at  Drexel  University 
[7-10]  to  study  the  effect  of  sequential  projectiles  Impacting  the  bottom 
of  a pre-drllled  hole  In  a target.  These  experiments  showed  In  general 
that  the  total  penetration  decreased  as  the  gap  between  the  projectiles  In- 
creased. "Post  mortem"  examination  of  the  targets  Indicated  that  the 
probable  cause  for  this  was  the  Interference  between  the  trailing  projec- 
tile and  the  occlusion  or  "back  flow"  caused  by  the  leading  projectile. 


In  this  report,  a numerical  study  of  sequential  Interior  impact  us- 
ing the  HELP  code  [11]  is  presented.  The  goal  of  this  study  is  to 
determine  the  mechanisms  involved  in  the  sequential  impact  process.  The 
experiments  of  [8-10],  which  use  low  speed  (1  mm/ps)  lead  projectiles  and 
lead  targets,  are  simulated  by  numerical  calculations.  The  numerical 
computation  of  copper  projectiles,  having  dimensions  of  typical  shaped 
charge  jet  segments,  impacting  at  typical  jet  velocities  onto  steel  tar- 
gets is  also  presented  in  this  study. 

In  Section  II,  the  results  of  the  breakup  study  are  presented.  The 
sequential  projectile  impact  calculations  are  given  in  Section  III. 
Finally,  conclusions  and  recommendations  for  future  study  on  these  topics 
are  presented  in  Section  IV. 
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II. 


SHAPED  CHARGE  JET  BREAKUP 


In  this  section,  the  details  of  our  current  study  on  shaped  charge 
jet  breakup  will  be  presented.  A new  one-dimensional  stability  theory 
will  be  derived  and  comparison  with  two-dimensional  numerical  calcula- 
tions will  be  made.  Additional  factors  affecting  jet  stability  using 
two-dimensional  calculations  will  also  be  presented. 

A.  Background 

In  a typical  shaped  charge  with  a conical  metallic  liner,  the  en- 
suing jet  travels  fastest  at  the  tip  with  the  velocity  decreasing  in  an 
approximately  linear  fashion  to  the  rear  of  the  jet  [12].  This  velocity 
gradient  causes  the  jet  to  stretch.  Because  of  this  stretching,  a terles 
of  surface  disturbances  or  "necking"  regions  eventually  develops  aicng  the 
length  of  the  jet.  Finally,  these  necked  regions  fail  and  the  jet  breaks 
into  a series  of  segments  [13-14].  In  this  report  we  will  study  the 
stability  and  growth  of  such  disturbances  in  shaped  charge  jets. 

In  the  study  of  flash  radiographs  of  metallic  jets,  we  have  observed 
that  there  exists  a remarkable  resemblance  between  the  breakup  process 
in  shaped  charge  jets  and  that  in  liquid  jets.  This  resemblance  is  re- 
ported in  [4]  and  [15].  It  is  a well  known  fact  that  for  liquid  jets  this 
breakup  is  caused  by  a surface  instability.  The  subject  of  liquid  jet 
instability  has  been  investigated  by  many  researchers  over  the  past  century 
[16-21],  These  investigations  have  shown  that  for  a given  liquid  jet 
diameter  there  exists  a critical  wavelength  of  the  surface  disturbance 
which  causes  the  amplitude  of  the  disturbance  to  become  unstable  and  grow 
without  bound.  The  primary  driving  force  for  this  instability  is  the 
surface  tension  of  the  liquid. 

There  are  some  very  significant  differences  between  the  classical 
liquid  jet  breakup  and  shaped  charge  jet  breakup,  however.  In  the  metallic 
shaped  charge  jets  surface  tension  is  not  an  important  factor.  The  major 
driving  force  for  the  disturbance  growth  is  caused  by  the  strength  of  the 
jetting  metal.  In  addition,  the  shaped  charge  jets  are  stretching  at  very 
high  rates  of  strain.  Most  of  the  classical  analyses  of  liquid  jets  deal 
with  non-stretching  jets.  The  studies  which  do  treat  stretching  [22,23] 
neglect  the  inertia  terms  in  the  governing  equations.  This  approximation 
may  be  valid  for  relatively  low  strain  rates,  such  as  those  encountered 
in  drawing  threads.  However  for  the  rates  encountered  in  shaped  charge 
jets,  the  effect  of  inertia  is  of  significant  importance.  Note  that  a 
direct  extension  of  the  stretching  theories  of  Tomotlka  [22]  or  Mikami  [23] 
to  include  inertia  terms  and  Jet  strength  is  a difficult  and  cumbersome 
task.  This  is  due  to  their  formulation  which  is  essentially  a two-dimen- 
sional axlsymmetrlc  spatial  (Eulerlan)  formulation  that  is  later  reduced 
to  one  dimension  by  assuming  a periodic  solution  along  the  jet  length. 

In  previous  reports  [3,4],  we  have  demonstrated  that  the  concepts  of 
surface  instability  may  be  applied  to  elastic-plastic  shaped  charge  jets 
through  the  use  of  the  two-dimensional  finite  difference  calculations. 


Various  factors  such  as  yield  strength,  density,  disturbance  wavelength, 
and  random  disturbances  were  considered.  The  major  results  of  the 
stability  study  presented  In  [4]  may  be  briefly  stated  as: 

1.  Jets  with  higher  yield  strengths  break  sooner,  all  else 
being  equal. 

2.  Jets  with  lower  densities  will  break  sooner,  all  else 
being  equal. 

3.  For  the  shaped  charge  jet  calculated,  there  is  a critical 
time  for  the  growth  of  the  disturbance  amplitude.  Distur- 
bances introduced  early  do  not  grow  appreciably  before 
this  tl.-ne,  but  grow  rapidly  after  this  critical  time.  Dis- 
turbances introduced  after  this  time  also  grow  rapidly. 

4.  A critical  wavelength  (or  a range  of  wavelengths)  exists; 
disturbances  having  this  wavelength  grow  faster  than  all 
others.  The  length  of  the  broken  jet  segment  caused  by 
this  crlcical  wavelength  is  in  the  range  of  measured  jet 
segment  lengths. 

5.  When  irregular,  or  random,  surface  disturbances  are  intro- 
duced, the  growth  is  again  dominated  by  the  disturbance 
component  with  the  critical  wavelength.  The  jet  surface 
grows  into  a shape  similar  to  that  obtained  if  only  the  wave 
with  the  critical  length  were  Introduced. 

Further,  in  a recent  paper  [15],  we  have  also  presented  a one-dimen- 
sional stability  theory  for  a plastic  stretching  jet.  This  theory  pre- 
dicts all  of  the  pertinent  stability  effects  indicated  by  the  two- 
dimensional  calculations  except  that  of  critical  wavelength. 

In  the  present  study,  two  approaches  will  be  applied  to  the  problem 
of  shaped  charge  jet  stability.  The  first  is  an  improved  one-dimen- 
sional theoretical  analysis  which  will  predict  a critical  wavelength. 

This  is  accomplished  by  incorporating  a stress  distribution  along  the 
axis  of  the  jet  which  accounts  for  the  "stress  concentration"  at  the 
"necked"  region.  The  derivation  of  this  Improved  theory  is  presented 
in  the  next  sub-section. 

The  second  approach  is  the  continued  use  of  two-dimensional  finite 
difference  calculations  to  study  additional  effects  such  as  a random 
distribution  of  strength  or  velocity  throughout  the  jet,  elasticity, 
etc.,  as  well  as  for  comparison  with  one-dimensional  results. 

B.  One-Dimensional  Stability  Theory  for  Stretching  Plastic  Jets 
1.  Conservation  Equations 

For  the  one-dimensional  description  of  the  jet  motion,  it  is 
more  convenient  to  use  the  Lagrangian  coordinate.  The  two  independent 
variables  involved  are  the  time  variable  and  a variable  attached  to  the 
original  axial  position  of  the  mass  particles  in  the  jet.  Following 
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Courant  and  Friedrichs  [24],  we  define  a material  coordinate  a,  which 
assigns  a position  to  each  material  particle  of  the  Jet  and  remains 
with  that  particle  as  follows 


a 


^ rx(a,t) 
^0%  •'x(0,t) 


pA  dx 


(1) 


where  P Is  the  density  of  the  Jet,  A Is  Its  cross-sectional  area,  Pq 
Is  the  Initial  density,  Aq  Is  the  Initial  cross-sectional  area  of 
the  undisturbed  Jet,  and  x Is  the  spatial  coordinate  along  the  axis  of 
the  Jet  fixed  to  the  laboratory  In  the  Eulerlan  sense.  Note  that  the 
area  A Is  a function  of  both  location  x and  time  t;  this  Is  different 
from  the  one-dlmenslonal  flow  within  fixed  rigid  walls.  From  this 
definition  we  can  Immediately  write  the  continuity  equation,  or  con- 
servation of  mass  equation,  for  an  arbitrary  element  of  Jet  as 


pA/PqAo 


or 


pAx 


a 


(2) 


where  subscripts  a and  x Indicate  differentiation.  Assuming  that  the 
normal  component  of  stress  In  the  axial  direction  a Is  the  only  non- 
zero atress  during  the  motion  of  the  Jet  and  that  this  stress  a Is  the 
average  of  the  axial  stress  over  any  cross-sectional  area,  we  then  per- 
form a momentum  balance  over  an  arbitrary  element  of  Jet  to  obtain 

pA  X “ oA  + Ao  (3) 

tt  X X 


or.  In  terms  of  the  material  coordinate  a as  follows 


^0%  \t 


(aA)^ 


(4) 


where  the  subscript  t Indicates  material  time  differentiation  (l.e. 
following  a material  particle).  Equations  (2)  and  (4)  are  the  conser- 
vation of  mass  and  momentum  equations  In  terms  of  the  Independent 
variables  a and  t.  The  specification  of  a constitutive  equation  will 
complete  the  system  of  governing  equations. 


2.  Definition  of  Strain 


In  the  study  of  one-dlmenslonal  stretching  problems  (e.g.  the 
tension  test),  two  types  of  strain  are  typically  used  (see,  e.g.  Malvern 
[20]).  The  engineering  or  conventional  strain  Is  defined  as 

de  - dl/lQ  (5) 

where  de  represents  the  conventional  strain  increment,  dl  represents  the 
change  In  length,  and  Is  a reference,  or  Initial,  length.  This 
expression  Is  usiially  Integrated  from  to  t to  obtain  e - (I-IIq)IIq. 
Note  that  in  our  present  development  this  engineering  strain  may  be 
equivalently  defined  as  e ■ (dx-da)/da. 
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The  second  strain  definition  usually  encountered  is  the  logarith- 
mic, true  or  natural  strain.  This  is  defined  by 

dt  • di/i  (6) 


where  dc  is  the  natural  strain  Increment.  The  obvious  difference  is 
that  in  this  definition  the  deformed,  or  current,  length  is  used  as  a 
basis.  Upon  Integration,  we  obtain  e "•  fcn(£/io)t  where  tQ  again 
the  original  length.  Using  the  conservation  of  mass,  we  can  express 
tlie  natural  strain  as 


in 


p_A„ 

0 0 

pA 


For  an  incompressible  J-*! jrmat ion  we  have 


(7) 


c - 


(8) 


Since  most  empirical  stress-strain  functions  are  usually  expressed 
in  terms  of  natural  strain,  we  will  employ  the  definition  of  Eq.  (6)  in 
this  paper.  However,  we  may  have  just  as  easily  used  engineering  strain 
to  obtain  the  same  final  results. 


3.  Constitutive  Equation 


For  a one-dimensional  plastic  formulation,  it  is  sufficient 
to  prescribe  a uniaxial  stress-strain  curve  as  a description  of  the 
material  properties.  Let  us  consider  a universal  stress-strain  curve 
given  by  the  function 

T - T(c)  (9) 

where  T is  the  effective  true  stress  and  e is  the  effective  true  or 
natural  strain  (see  Malvern  [25],  pg.  364ff.  for  definitions  of  effec- 
tive stress  and  strain  to  be  used  with  the  Mlses  yield  condition).  In 
a truly  one-dimensional  stress  problem,  e.g.  uniaxial  tension,  the 
effective  stress  and  strain  simply  reduce  to  the  axial  components  of 
stress  and  strain  respectively.  Thus,  Eq.  (9)  may  be  obtained  from  a 
uniaxial  tension  test.  Some  older  authors  (see,  e.g.,  Bridgman  [26]) 
use  the  term  "flow  stress"  for  the  true  stress  when  the  other  stresses 
vanish.  Thus  the  stress-strain  law  given  by  Eq.  (9)  is  the  flow  stress 
as  a function  of  axial  strain  in  a one-dimensional  sense  when  no  neck- 
ing occurs.  At  present  we  will  not  specify  the  actual  form  of  Eq.  v.9). 

The  material  is  also  assumed  to  be  incompressible  ’'f  that 

p - pQ  ■ const.  (10) 

4.  Stretching  Jet  without  Surface  Disturbance 

Let  us  first  formulate  the  initial  and  boundary  conditions 
that  most  closely  describe  the  stretching  jet  of  the  shaped  charge  without 
surface  disturbances.  The  solution  of  this  problem  will  provide  a "base 
solution"  upon  which  we  can  superimpose  a "disturbance  solution"  to  study 
the  stability. 
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In  the  Jet  without  surface  disturbances,  It  Is  valid  to  assume 
that  the  stress  field  Is  truly  one-dlmens lonal  with  o - T(c).  Using 
this  fact  together  with  Eqs.  (A), (8), (9),  and  (10)  we  have  the  follow- 
ing final  form  of  the  momentum  equation 


<^0%  *tt 


(T  - ^)A  . 


(11) 


Equations  (2)  and  (11)  are  the  governing  equations  for  the  one-dlmen- 
slonal  stretching  problem  without  surface  disturbances  In  terms  of  the 
dependent  variables  x and  A. 


We  may  assume  that  one  end  of  the  jet  Is  fixed  without  any  loss 
of  generality.  The  other  end  of  the  Jet  Is  assumed  to  move  at  a con- 
stant velocity  Vq.  Initially,  the  velocity  varies  linearly  between  the 
ends.  The  Initial  cross  sectional  area  of  the  Jet  Is  Aq,  a constant; 
the  Initial  length  of  the  Jet  Is  Lq.  It  Is  convenient  to  define  an 
Initial  stretching  rate  as  hq  “ ^q/Lq.  This  configuration  is  depicted 
In  Fig.  1.  Mathematically  we  express  these  boundary  conditions  as 


Xt  - 0 

at 

a - 0 

(12) 

at 

a - Lq 

(13) 

and  the  Initial  conditions  as 

Xt  - noa 

at 

t - 0 

(lA) 

A - Aq 

at 

t ="  0. 

(15) 

Tc  can  be  shown  that  the  solution 

X - a(nQt 

+ 1) 

(16) 

A - A-/(n, 

,t  + 1) 

(17) 

satisfies  the  governing  equations,  Eqs.  (2)  and  (11),  and  the  auxiliary 
conditions  Eqs.  (12)-(15).  Equations  (16)  and  (17)  represent  the  solu- 
tion to  the  stretching  Jet  problem  without  surface  disturbances.  Note 
that  this  base  solution  Is  Independent  of  the  constitutive  equation  and 
thus  Eqs.  (16)  and  (17)  are  valid  for  any  Incompressible  matarlal.  Also, 
all  Jet  particles  maintain  constant  velocity  at  their  Initial  value, 
with  no  acceleration,  or, 

x^  “ riQfl  for  t ^ 0.  (18) 

and  the  Jet  area  decreases  with  time  but  Is  uniform  at  any  given  time. 

5.  Stretching  Jet  with  Initial  Disturbances 

Following  the  familiar  approach  In  studying  the  hydrodynamic 
stability  of  Jets,  we  will  Introduce  a small  disturbance  upon  the  base 
flow  and  find  the  solution  to  the  problem  with  this  disturbed  initial 
condition.  If  the  disturbance  dies  down  as  time  progresses,  the  flow 
Is  stable;  if  It  grows  with  time,  the  flow  Is  unstable.  Because  of  the 
difficulty  Involved  in  treating  nonlinear  equations,  we  will  consider  the 
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linearized  problem,  which  Is  valid  for  disturbances  of  small  amplitude, 
and  for  small  strains  after  disturbance  Initiation. 


Once  a surface  disturbance  Is  Introduced  In  the  Jet,  the  stress 
solution  Is  no  longer  exactly  one-dlmenslonal,  but.  In  fact,  the  average 
axial  stress  o varies  along  the  axis  of  the  Jet.  For  a static  tensile 
test,  Bridgman  [26]  developed  the  three-dimensional  solution  In  the  region 
of  the  neck.  His  result  showed  that  the  axial  stress  at  the  neck  was 
equal  to  the  flow  stress  at  the  free  surface  and  Increased  to  a maximum 
on  the  axis.  Thus,  the  average  axial  stress  at  the  neck  is  actually 
higher  than  the  flow  stress.  Bridgman  obtained  his  results  by  using  the 
Mlses  yield  criterion  and  the  assumption  of  Isotropic  hardening.  This 
Is  equivalent  to  postulating  a universal  stress-strain  law  such  as  Eq. 

(9)  (see  Malvern  [25],  pg  364).  The  result  obtained  by  Bridgman  for  the 
average  axial  stress  Is  given  by 

0 . T^i  + 2f)  in  (1  + i f)  (19) 

n 

where  r Is  the  outside  radius  at  the  neck  and  R Is  the  radius  of  curva- 
ture of  the  neck. 

We  now  postulate  the  use  of  the  Bridgman  stress  solution  for  our 
stretching  Jet  problem  at  any  cross-section  (l.e.  not  only  at  the  neck). 
Further  we  restrict  our  result  to  small  disturbances  which  Implies  that 
1/R  S 3^r/3x^.  Under  these  conditions  we  have  for  the  average  axial 
stress  distribution  the  following 

2 

0 - T(1  + J (20) 

3x 

where  r • r(x,t)  Is  the  radius  of  outer  disturbed  Jet  surface.  Com- 
blnl.ng  Eqs.  (4) , (8) , (9) , (10) , and  (20)  we  obtain  the  following  form  of 
the  momentum  equation 

Note,  in  writing  the  final  form  of  this  equation  we  have  dropped  higher 
order  terms  in  the  derivatives  of  A (e.g.  A^,  A A , etc.).  This  is 
valid  since  we  eventually  will  restrict  our  solution  to  small  surface 
disturbances  and  as  will  be  seen  these  terms  will  only  involve  higher 
order  disturbance  quantities.  Also  note  that  If  no  surface  disturbance 
exists  the  second  term  on  the  right  hand  side  of  Eq.  (21)  vanishes  and 
Eq.  (21)  simply  reduces  to  Eq.  (11). 

The  governing  differential  equations  are  Eqs.  (2)  and  (21),  and 
the  boundary  conditions  Eqs.  (12)  and  (13)  will  remain  the  same.  We 
will  first  specify  the  disturbed  initial  condition  for  the  area  A only, 
and  will  not  Impose  any  specific  Initial  condition  for  x. 
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Consider  the  following  disturbed  Ini'^lal  condition  for  the  area, 
A(a,0)  = ^0^^  ^ ^0 


where 


b 


2irn 


n Is  the  number  of  waves  along  the  given  jet  length,  Xq  Is  the  Initial 
wavelength  of  the  disturbance,  and  Iq  Is  the  Initial  amplitude  of  the 
disturbance  function  which  Is  prescribed  to  be  a small  quantity.  It  will 
be  shown  that  the  following  area  function  can  satisfy  the  governing 
equations  and  has  the  Initial  value  given  by  Eq.  (22) 

A(a,t)  - Aq[1  + f(t)co8  baJ/Cript  1)  (23) 


where  f(t)  Is  a function  governing  the  amplitude  of  the  disturbance 
term  and  f (0)  “ Iq.  Substituting  this  expression  Into  Eq.  (2),  we 
obtain 


X - (n«t  + 1)/[1  -*■  f(t)co8  ba].  (24) 

a u 

If  we  neglect  terms  of  the  order  f^  and  higher,  then  Eq.  (24)  may  be 
written  as 


X ■■  (n/>t  + 1)[1  - f(t)cos  ba].  (2S) 

a u 

Integrating  this  expression  and  using  the  boundary  condition,  we  obtain 
for  the  position 


X “ (HqI  + l)[a  - (f/b)sln  ba]. 


(26) 


Substituting  the  expressions  for  x and  A,  Eqs.  (26)  and  (23),  and  their 
derivatives  Into  Eq.  (21)  and  neglecting  terms  of  the  order  f^  and 
higher,  we  obtain 


ftt  * 


V 


4ir  PQ(nQt+l)' 


f-0- 


(27) 


If  the  function  f satisfies  Eq.  (27),  then  the  corresponding  Eqs.  (23) 
and  (26)  are  the  solution  of  the  linearized  problem  of  Eqs.  (2)  and 
(21),  botindary  conditions  Eqs.  (12)  and  (13),  and  the  disturbance  Initial 
condition  Eq.  (22).  The  other  initial  condition  on  x^  turns  out  to  be 

x^  - HqS  - ^lhQf(O)  + f^(0)]sln  ba  at  t - 0.  (28) 


Up  to  this  point  we  have  not  considered  the  form  of  the  stress- 
strain  law  Eq.  (9).  We  will  now  consider  two  cases,  simple  work  harden- 
ing and  perfectly  plastic. 
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Case  I - Work  Hardening 


With  no  loss  of  generality,  we  will  consider  the  strain  hardening 
law  to  be  a series  expansion  about  the  initial  point,  e ■ 0,  T ■ o^, 

T(e)  - Oq  + Ce  + Fe^  + •••  (29) 

Since  we  seek  a linearized  solution,  we  may  restrict  our  result  to  small 
strains  from  the  disturbance  initiation.  Under  this  small  strain 
assumption. 


T 


r*  a. 


il 

dc 


C . 


(30) 


We  further  neglect  nQt  with  respect  to  1 in  the  coefficient  of  f in 
Eq.  (27)  which  is  consistent  with  the  small  strain  assumption*.  These 
approximations  applied  to  Eq.  (27)  yield  the  following  differential 
equation  governing  the  disturbance  function 


2no(not+l)f^  - — 


2 a 


-4 


0 (31) 


Equation  (31)  is  an  ordinary  differential  equation  of  the  Euler 
type  and  can  be  solved  by  elementary  methods.  The  solution  is  of  the 
form 


f(t)  - Q(nQt  + D' 

where  Q is  an  arbitrary  constant,  and 


8 “ -I  1 ± •'I  + 

■ 4 fe)  ■ CV)  ] ■ 


(32) 


(33) 

(34) 


The  disturbance  function  f(t)  will  decay  with  time  if  s is 
negative  and  grow  if  s is  positive.  From  Eqs.  (33)  and  (34)  it  can 
be  seen  that  a positive  root  of  s exists  if 


(35) 


* Note  that  we  do  not  neglect  no^  in  the  other  two  terms  since  it  is  not 
necessary  to  obtain  a simple  solution  and  no  problem  is  encountered  by 
carrying  it. 
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In  our  analysis  so  far,  we  have  assumed  that  there  is  an  initial  state 
where  the  strain  Is  zero  and  the  length  of  the  jet  is  Lq.  The  stress- 
strain  relation  is  expanded  into  a series  about  this  state.  This  initial 
state,  or  Initial  time,  is  completely  arbitrary.  We  may  consider  any 
t Ime  during  the  Jet  stretching  process  as  the  initial  time  and  the  sta- 
bility analysis  presented  above  is  always  true.  In  this  sense,  oq  in 
Eq.  (29)  is  the  flow  stress  at  the  time  under  consideration  and  C is  the 
slope  (or  tangent  modulus)  of  the  stress-strain  curve  at  that  stress  level  . 
The  stability  results,  Eqs.  (32),  (33),  and  (34)  depend  only  on  the  values 
of  the  stress  and  the  slope  of  stress-strain  curve,  and  are  Independent 
of  the  value  of  the  strain.  Note  also  that  in  our  derivation  the  stress 
is  the  true  stress,  or  force  per  unit  current  area,  and  the  strain  is  the 
natural  or  true  strain.  In  contrast  to  the  more  elementary  theory  [15] 
which  did  not  account  for  the  stress  concentration,  the  stability  results 
do  depend  on  the  initial  dimensions  tq  and  Xq.  We  thus  see  that  there 
is  an  initial  wavelength  which  yields  a maximum  disturbance  growth.  We 
define  this  to  be  the  critical  wavelength.  By  differentiating  Eq.  (33) 
with  respect  to  Xq,  holding  all  other  parameters  fixed,  and  setting  the 
result  to  zero  we  find  the  critical  wavelength  to  be  given  by 


(36) 


We  further  note  that  for  a critical  disturbance  given  by  Eq.  (36),  the 
condition  for  a positive  exponent  Eq.  (35)  reduces  simply  to  Oq  > C, 
l.e.  if  the  stress  is  greater  than  the  tangent  modulus. 


We  may  now  state  the  jet  stability  criterion  as  follows: 

The  stability  of  a stretching  jet  may  vary  from  time  to  time 
For  a jet  having  a critical  disturbance  defined  by  Eq.  (36), 
when  the  true  stress  Is  less  than  or  equal  to  the  slope  of 
the  true  stress-true-straln  curve,  the  jet  Is  stable;  other- 
wise the  jet  is  unstable.  That  is,  a jet  having  a critical 
disturbance  is  stable  if 


Oq  1 C.  (37) 

The  stability  of  a jet  having  a disturbance  other  than  critical 
is  governed  by  the  criterion  of  Eq.  (35).  For  an  unstable  jet, 
any  small  disturbance  will  grow,  and  the  rate  of  growth  is 
governed  by  the  value  of  K. 


It  is  interesting  to  note  chat  in  Che  conventional  static  simple 
tension  tests,  the  plastic  Instability  criterion  [27]  is  Che  same  as 
our  present  one  for  critical  disturbances  although  completely  different 
physical  reasoning  was  used.  In  the  static  case,  the  critical  insta- 
bility point  la  defined  as  the  point  where  the  engineering  stress  vs. 
engineering  strain  curve  has  a zero  slope.  This  point,  when  transformed 
on  the  true  stress  vs.  true  strain  curve,  is  the  point  where  ■ C. 
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The  argument  for  establishing  this  static  criterion  Is  that  when  the 
engineering  stress  Is  a maximum,  the  load  carrying  capacity  of  the 
specimen  reaches  a maximum.  The  specimen  will  stretch  to  complete  fail- 
ure while  the  total  load  Is  decreasing. 

The  present  disturbance  solution  yields  some  very  Interesting 
trends.  From  Eqs.  (33)  and  (3A)  it  can  be  seen  that  the  larger  the 
value  of  K,  the  faster  the  growth  rate.  This  factor  is  strongly  depen- 
dent upon  the  ratio  Oq/pq.  Thus,  we  see  that,  in  general,  for  unstable 
jets  the  larger  the  flow  stress,  the  faster  the  growth  rate;  the  larger 
the  density,  the  smaller  the  growth  rate.  It  is  Interesting  to  note 
that  the  growth  rate  is  inversely  proportional  to  the  initial  stretch- 
ing rate  Hq;  as  the  stretching  rate  becomes  larger  the  growth  becomes 
slower. 

Case  II  - Perfectly  Plastic 


For  a perfectly  plastic  material,  Eq.  (9)  becomes  independent  of 
strain  and  is  given  by 


T(e)  - Y 


(38) 


where  Y is  the  yield  stress  in  simple  tension  and  taken  as  a constant 
for  a given  material.  Using  Eq.  (38)  together  with  Eq.  (27)  and  again 
dropping  riQt  with  respect  to  1 in  the  f term  coefficient,  we  obtain  a 
differential  equation  similar  to  Eq.  (31).  The  solution  is  again  given 
by  Eqs.  (32)  and  (33)  but  now  the  definition  of  K,  Eq.  (34),  is  given  by 


(39) 


We  again  note  the  existence  of  a stability  criterion  based  on  whether 
the  exponent  s is  positive  or  negative.  However  for  perfectly  plastic 
Jets,  the  stability  is  not  based  on  material  properties  but  only  on  the 
initial  geometry  of  the  disturbance.  From  Eq.  (39)  we  see  that  the  jet 
is  stable  if 


< ii/2 


otherwise  the  jet  is  unstable. 


(40) 


Also,  from  Eqs. 
critical  wavelength, 
ing  case  we  find  the 
be  given  by 


(32),  (33)  and  (39)  we  note  the  existence  of  a 
Following  the  procedure  used  in  the  work  harden- 
crltlcal  wavelength  for  perfectly  plastic  jets  to 


(41) 


Thus,  in  the  perfectly  plastic  case  the  initial  critical  wavelength  is 
always  2.22  times  the  initial  Jet  diameter. 
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Again  the  same  growth  trends  appear  In  the  perfectly  plastic  case 
as  was  observed  in  the  work  hardening  case. 

6.  Complete  Stability  Solution 

In  Section  IID,  we  will  compare  this  linearized  theory  for  the 
case  of  a perfectly  plastic  material  with  a two-dimensional  calculation 
for  a particular  problem  with  specific  initial  conditions.  For  this 
purpose,  a complete  solution  to  Eq.  (31)  is  needed.  The  general  solu- 
tion of  Eq.  (31)  is 

S 8 

f(t)  - Qi(nQt  + 1)  ^ + Q2(not  + D ^ (^2) 

where  Qj^  and  are  constants  to  be  determined  from  Initial  conditions 
on  f,  and  sj  and  S2  are  the  two  roots  defined  by  Eq.  (33).  For  this 
specific  problem  we  Impose  initial  conditions  in  which  the  position  of 
the  surface  is  disturbed  and  at  rest.  These  can  be  approximated  as 

f(0)  - f^ 
f^(0)  - 0. 

Using  these  two  conditions,  we  obtain  the  following  result 

f S3 

“ 7^  ^^2%^  + 1)  ^ - s^(not  + 1)  2] 

This  disturbance  function  can  now  be  used  to  compute  the  position, 
velocity,  and  cross-sectional  area  for  the  stretching  jet  problem  with 
a surface  disturbance. 

C.  Two-Dimensional  Finite-Difference  Calculation 
of  Stretching  Elastic-Plastic  Jets 

The  numerical  study  of  jet  stability  was  undertaken  through  the 
use  of  the  two-dimensional  code  HEMP  [5].  The  HEMP  code  is  a general 
purpose  code  which  solves  the  conservation  equations  of  two-dimensional 
elastic-plastic  flow  in  plane  coordinates  or  in  axlsytnmetrlc  coordinates. 

The  solution  is  by  the  method  of  finite  differences  and  uses  the 
Lagrangian  formulation.  The  code  has  the  capability  of  handling  many 
various  boundary  and  initial  conditions.  In  this  section  we  will  briefly  < 

describe  the  general  approach  of  applying  HEMP  to  the  jet  stability  prob- 
lem. 

Karpp  [6]  first  applied  HEMP  code  calculations  to  the  problem  of 
a stretching  elastic-plastic  jet  with  the  surface  slightly  disturbed. 

After  calculating  surface  disturbances  of  various  wavelengths,  he  found  f 

a range  of  most  unstable  wavelengths  with  reasonable  growth  rates.  We  i 

have  applied  the  same  basic  method  to  study  the  breakup  problem  in  more  f 

detail.  We  have  used  the  HEMP  code  to  study  many  various  aspects  of 
stretching  elastic-plastic  Jets  including  the  effects  of  yield  strength,  , 

disturbance  initiation  time,  disturbance  wavelength,  random  surface  dis-  i 

turbances,  and  inertia  forces. 
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(43) 


(44) 
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For  these  two-dimensional  calculations,  a stretching  shaped  charge 
]et  Is  modelled  by  an  axisymmetric  bar  fixed  at  one  end,  with  the  other 
end  moving  at  a constant  velocity.  A linear  velocity  distribution  in 
the  axial  direction  is  again  imposed  as  the  initial  condition,  and  the 
surface  of  the  bar  is  initially  perturbed  in  the  shape  of  a cosine 
function,  as  previously  shown  in  Fig.  2.  Note  that  to  a first  order 
approximation  perturbing  the  radius  is  equivalent  to  perturbing  the  cross- 
sectional  area  as  done  in  the  one-dimensional  theory  Eq.  (22).  The 
perturbed  surface  is  free  from  any  tractions  and  the  end  surfaces  are 
free  of  any  shear  stresses.  Let  the  axial  velocity  be  v(x,r,t),  the 
radial  velocity  be  u(x,r,t),  the  stress  vector  on  the  lateral  surface  be 
o(x,rg,t),  and  the  shear  stresses  on  the  end  surfaces  be  T^^(0,r,t), 

T^Q(0,r,t),  T^^(LQ,r,t),  T^Q(LQ,r,t),  then  the  boundary  conditions  are 

v(0,r,t)  • 0 
v(L,r,t)  - Vq 
j2(x,r^,t)  - Q 

Trx(O.r.t)  - T^Q(0,r,t)  - T^^(LQ,r,t)  - T^Q(LQ,r,t)  - 0 

and  the  initial  conditions  are 


v(x,r,0) 


(Vp/Lq)x  and  u(x,r,0) 


where  r^  is  the  disturbed  surface  tq  + D-(cos  2Tinx/LQ) . The  calculation 
of  one  cycle  of  the  surface  wave  along  the  axial  direction  is  sufficient 
because  of  the  symmetry  of  the  problem. 

Using  this  basic  approach  we  have  calculated  various  problems  of 
stretching  jets.  Unless  otherwise  noted,  we  have  considered  elastic- 
perfectly  plastic  copper  jets  with  an  equation  of  state  for  copper  as 
specified  in  the  HEMP  code.  The  value  of  density  used  was  8.9  Mg/m^; 
the  elastic  shear  modulus  was  taken  as  45.6  GPa  and  the  yield  stress 
used  was  0.2  GPa. 

D.  Comparison  of  One-D  and  Two-D  Results 

To  compare  the  results  of  the  linearized  one-dimensional  theory  and 
the  two-dimensional  numerical  solutions,  we  have  computed  several  prob- 
lems having  dimensions  and  stretching  rates  which  span  the  range  of  those 
normally  encountered  in  typical  shaped  charge  jets.  More  specifically, 
we  have  considered  two  stretching  rates  in  these  problems  no  “ 0.029  us"^ 

and  Hq  ~ 0.008  ps  In  all  cases  the  initial  undisturbed  radius  was 
rQ  ■ 1.5mm.  All  jets  were  considered  to  be  perfectly  plastic  with  mate- 
rial constants  specified  in  Section  IIC.  The  wavelength  Xq  was  varied 
such  that  1.5  < Xq/ZTq  < 6.  In  our  comparisons  we  will  use  the  relative 
growth  A defined  as 


A(t)  - lD(t) 
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where  0(1)  Is  the  amplitude  of  the  surface  disturbance  at  time  t and  Dq 
Is  the  Initial  amplitude.  In  terms  of  the  disturbance  function  f,  it 
can  be  shown  that  for  small  disturbances  (l.e.  small  f) 


Substituting  the  one-dlmenslonal  solution  Eq.  (44)  into  Eq.  (48),  we 
see  that  Iq  Is  eliminated  and  thus  A for  small  disturbances  is  inde- 
pendent of  the  Initial  disturbance  Eq.  Figure  2 shows  a plot  of 
A vs  XqHxq  at  a time  of  28  us  after  disturbance  initiation  for  both 
strain  rates  considered.  In  this  case  the  two-dimensional  solution  had 
Dq  « 0.001  Cq.  The  two-dimensional  solutions  show  that  the  critical 
length  to  diameter  ratio  (l.e.  the  one  which  yields  maximum  A)  varies 
with  the  initial  strain  rate  from  1.5  at  Hq  = 0,029  us~^  to  2.2  at 
Pq  “ 0.008  us”^.  As  discussed  above,  for  perfectly  plastic  jets,  the 
one-dimensional  theory  always  predicts  crit  * regardless  of 

strain  rate.  We  should  also  note  that  Karpp  and  Simon  [6]  found 

“ 2.2  using  HEMP  code  for  jets  with  strain  rates  in  the 
0 0 crlt  ® 

range  of  0.005  us  Dq  < 0.008  us  Further,  the  growth  predicted  by 
the  one-dimensional  theory  appears  in  good  agreement  with  the  two- 
dimensional  solution  at  the  lower  strain  rate.  At  the  higher  strain 
rate,  0.029  us“^,  however,  the  theory  predicts  values  which  are  quite 
low.  To  exemplify  this  further  Figs.  3 and  4 show  plots  of  A vs  t for 
Pq  “ 0.029  us”^  and  pq  - 0.008  us“^  respectively.  Both  have  an  initial 
length  to  diameter  ratio  of  2.2.  Also  in  Fig.  3 we  show  a two-dimen- 
sional solution  with  the  initial  amplitude  increased  by  a factor  of  50 
to  Dq  « .05  tq.  We  see  that  for  this  set  of  conditions  the  relative 
growth  is  also  independent  of  initial  amplitude.  It  is  also  important 
to  note  that,  for  a given  time  after  disturbance  initiation,  the  linear- 
ized solution  for  the  case  with  higher  initial  strain  rate  involves  more 
error  than  for  the  case  with  lower  initial  strain  rate.  It  is  antici- 
pated that  a solution  which  uses  the  linearized  solution  for  successive 
time  steps  will  alleviate  this  problem. 


It  is  interesting  to  note  that,  for  the  small  amplitude  problem, 

Dq  - 0.001  rQ,  the  HEMP  code  solution  Indicates  that  the  material  is 
completely  in  the  plastic  range.  After  some  initial  fluctuation  in  the 
solution,  HEMP  predicts  an  axial  stress  at  the  prescribed  yield  point 
0.2  GPa  with  the  other  stresses  being  negligible.  Therefore,  neglecting 
elastic  unloading  in  the  theory  is  a reasonable  assumption.  For  problems 
with  larger  disturbance  amplitudes,  or  if  we  continue  the  present  calcu- 
lation until  the  disturbance  grows  significantly  larger,  the  HEMP  solu- 
tion does  not  remain  completely  plastic  and  an  elastic  region  forms  near 
the  ends,  where  the  cross-sectional  area  is  larger.  All  of  the  stretch- 
ing then  takes  place  in  the  "necked"  region.  We  will  discuss  the 
effects  of  disturbance  amplitude  and  elasticity  in  more  detail  in  the 
next  section. 
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Figure  2.  Comparison  of  relative  growth  predicted  by  one-dimensional  theory  and 

two-dimensional  calculations  for  two  different  strain  rates  and  various 
length  to  diameter  ratios  at  28  ps  after  disturbance  initiation. 
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We  have  noted  above  that  for  a perfectly  plastic  jet  the  theoret- 
ical solution  depends  on  the  ratio  Y/Pq-  To  observe  If  this  result  Is 
In  agreement  wltl>  HEMP  predictions,  we  have  computed  a series  of  problems 
which  were  Identical  except  the  Individual  values  of  Y and  Pq  were 
varied  keeping  the  ratio  Y/Pq  constant.  We  varied  Y and  p from  1/3  to 
2 times  the  values  used  for  copper  in  the  present  problem,  and  the  solu- 
tion for  the  minimum  cross-sectional  area  did  not  deviate  more  than  1%. 
Also,  we  have  found  that  as  the  ratio  Y/pq  Increases,  the  disturbance 
growth  increases  (see,  e.g,  [4))  which  Is  also  In  agreement  with  the 
theoretical  trends. 

Further,  we  would  like  to  note  that  the  trends  predicted  by  the 
theory  are  also  In  general  agreement  with  the  numerical  calculations 
presented  by  Karpp  and  Simon  [6],  They  have  also  Indicated  that  growth 
increases  with  Increasing  yield  strength  and  decreases  with  increasing 
stretching  rate. 

Finally,  we  have  performed  numerical  computations  with  a work 
hardening  version  of  HEMP  recently  developed  by  Karpp  [28].  These  cal- 
culations show  that  the  presence  of  a small  amount  of  hardening  during 
the  flow  causes  the  disturbance  to  stabilize  and  decay  in  agreement 
with  the  stability  criterion  set  forth  in  Section  IIB.  Details  of  these 
calculations  will  be  reported  in  the  future. 

In  the  next  sub-section,  we  will  further  explore  the  use  of  two- 
dimensional  calculations  to  study  various  other  effects  on  jet  stability. 

E.  Additional  Two-Dimensional  Stability  Results 

Two-dimensional  calculations  of  the  stability  of  stretching  jets 
were  used  extensively  during  the  course  of  this  study.  In  addition  to 
the  two-dimensional  results  presented  above,  various  other  effects  were 
studied  using  the  HEMP  code.  The  results  of  these  calculations  will  be 
presented  in  this  sub-section. 

1.  Estimation  of  Jet  Yield  Strength  Using  Measured 
Breakup  Time 

From  the  measurement  of  broken  jet  radiographs  (see  [1]  and 
[4])  we  have  determined  that,  for  many  jets  (e.g.  the  jet  from  the 
81.3mm  BRL  standard  charge  [1]),  the  breakup  time  varies  along  the 
length  of  the  jet.  If  we  hypothesize  that  the  yield  strength  varies  along 
the  jet  and  causes  the  variation  in  breakup  time,  we  may  then  determine 
this  variation  in  yield  strength  by  HEMP  code  calculations,  if  the 
measured  breakup  time  data  are  used.  To  do  this,  we  have  selected  three 
elements  along  a typical  jet  (81.3mm  BRL.  Standard  charge  [1])  which  have 
three  different  breakup  times.  Let  us  denote  these  elements  as  A,  B, 
and  C.  These  elements  will  break  into  segments  at  times  of  t/^)  ~ 94  us, 

(B)  (C)  ^ 

t^  - 121  us  and  - 161.6  us  as  determined  from  the  least  squares 

breakup  distribution  of  [1].  These  times  are  based  on  t ~ 0 when  the 
detonation  wave  reaches  the  cone  apex.  These  elements  have  average  jet 
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4.9  km/s. 


velocities  of  - 7.1  km/s,  « 6.2  km/s,  and  ■ 

Using  a one-dimensional  model  we  may  trace  back  these  elements  to  their 
original  position  on  the  cone  and  also  compute  their  impinging  velocity 
onto  the  stagnation  point  during  formation.  From  these  one-dimensional 

calculations  we  find  ” 0.44,  ■*  0.51,  and  x^^^  » 0.60  and 

- 0.85,  » .75,  and  = 0.62,  where  x is  the  axial  position 

from  the  cone  apex  divided  by  the  cone  height  (h)  and  M is  the  Mach 
Number  of  the  impinging  flow. 

Now,  we  consider  each  of  these  elements  to  be  stretching  before 
breakup  and  perform  HEMP  calculations  of  each.  We  introduce  a sinusoidal 
surface  disturbance  as  described  in  Section  IIC  in  each  element  at  a 
considerable  time  before  the  breakup  time  and  observe  if  the  disturbance 
growth  will  cause  the  element  to  break  at  the  desired  breakup  time.  A 
schematic  position-time  plot  of  this  process  is  shown  in  Fig.  5.  The 
yeild  strength  in  each  element  is  adjusted  until  the  desired  result  Is 
obtained.  The  HEMP  code  does  not  have  the  capability  to  actually  pre- 
dict a break  in  the  stretching  element.  We  prescribe  a somewhat  arbi- 
trary criterion  for  a "necked"  section  to  break  which  is  when  the  necked 
cross-sectional  area  is  10%  of  the  average  maximum  cross-sectional  area 
of  the  Jet  segments  as  measured  from  Jet  radiograph.  The  initial  con- 
figuration for  each  element  of  Jet  used  in  the  HEMP  calculations  is  given 
in  Table  I as  is  the  time  of  disturbance  initiation  tj^.  The  difference 
in  velocity  between  the  two  ends  of  each  element  A Vj  is  also  given  and 
was  computed  using  the  velocity  gradient  in  the  Jet  at  the  time  t^.  The 
results  of  the  HEMP  computations  are  shown  in  Table  II,  where  we  see 

that  strengths  of  - 0.12  GPa  and  - 0.066  GPa  and  Y^^^  = 0.036 

GPa  cause  breakup  times  very  close  to  the  experimental  values. 


Table  I 

Initial  configuration  for  segments  A,B,  and  C 
used  in  the  HEMP  calculations. 


' element 

property  __ 

A 

B 

C 

initial  length 
(mm) 

4.79 

4.65 

4.59 

initial  radius 
(mm) 

1.60 

1.55 

1.53 

(mm/us) 

.147 

.105 

.078 

i "i 

' (us) 

43.6 

54.8 

74.3 

initial  distur- 
bance amplitude  (mm) 

.080 

.078 

.077 
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Table  II 


Results  of  HEMP  calculations  and 
comparison  with  experimental  breakup  time. 


Element 

No. 

x/h 

M 

Exp. 

b.u.  time 

(ms) 

HEMP 

- 

A 

.44 

00 

94 

94 

0.12 

B 

.51 

.75 

121 

119 

0.066 

C 

.60 

.62 

161 

159 

0.036 

2.  Random  Yield  Strength 

We  have  examined  the  effects  of  surface  disturbance  through- 
out our  work  so  far.  Suppose  now  we  have  no  surface  disturbance,  but 
the  yield  strength  of  the  material  varies  from  cell  to  cell  In  the 
HEMP  calculations.  As  before,  we  have  used  a portion  of  jet  subject 
to  a linear  velocity  gradient,  which  Is  shown  In  Fig.  6.  Note  that 
only  the  portion  above  the  centerline  is  shown  for  convenience.  The 
calculation  grid  Is  divided  Into  96  cells  and  the  yield  strengths  of 
each  cell  are  distributed  randomly  throughout  the  grid,  assuming  the 
following  distribution: 


of  cells 

Yq  (GPa) 

40 

0.200 

15 

0.210 

15 

0.190 

10 

0.175 

10 

0.225 

3 

0.160 

3 

0.240 

The  initial  dimensions  and  velocity  gradient  are  chosen  such  that  the 
portion  of  jet  shown  will  be  equal  to  two  average  segment  lengths  of 
jet  at  the  breakup  time,  as  measured  from  the  radiographs.  Figure  6 
shows  the  configuration  of  the  jet  at  18  ws  and  36  ps.  It  is  Inter- 
esting to  note  that  It  appears  to  break  Into  two  full  segments  plus  a 
little  extra.  Thus  we  see  that  the  HEMP  calculation  predicts  segment 
length  In  the  neighborhood  of  the  measured  lengths. 

3.  Details  of  the  Stress  Distribution  In  a Jet 


In  the  Appendix  to  this  report  we  have  compiled  for  easy 
reference  the  stress  distribution  as  predicted  by  HE^  for  a stretch- 
ing jet  with  a sinusoidal  surface  disturbance.  The  dimensions  of  the 
initial  configuration  are  also  given  there.  Tabular  results  are  given 
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at  several  times  of  Interest  and  several  plots  of  axial  stress  versus 
time  are  Included.  Note  that  the  tabular  data  are  given  at  very  early 
times  and  at  several  times  Itv  the  neighborhood  of  when  the  first  appear- 
ance of  elastic  cells  occurs.  This  was  done  to  see  the  "start  up" 
stress  distribution  as  computed  by  HE31P.  The  later  times  were  tabulated 
to  observe  the  effects  of  the  existence  of  elastic  cells  on  the  growth 
of  the  disturbance. 

Finally,  the  stress  time  plots  presented  show  that  at  early  times 
the  axial  stress  Is  somewhat  erratic.  In  the  necked  region,  we  see  that 
the  cells  remain  plastic  and  the  stress  "levels  off"  at  approximately 
the  yield  stress  In  simple  tension,  which  Indicates  In  this  region  the 
flow  may  be  approximately  one-dlmenslonal.  In  the  region  away  from  the 
necked  region  we  see  that  the  cells  eventually  become  elastic  with  a 
general  decrease  In  the  axial  stress. 

4.  The  Effect  of  Calculatlonal  Mesh  Size 

We  have  examined  the  effect  of  using  fewer  grid  cells  In  the 
HEMP  code  calculations  of  a jet  stability  problem  (since  fewer  cells 
will  be  more  economical  to  calculate).  In  our  typical  calculations  we 
have  been  using  5 cells  across  the  radius  of  the  jet  and  16  cells  along 
each  disturbance  wavelength.  To  observe  the  effect  of  mesh  size,  we 
have  considered  two  problems  which  are  Identical  except  for  the  cell 
structure.  One  has  the  standard  5x16  structure  whereas  the  other  only 
has  2 cells  across  the  radius  and  6 cells  along  one  disturbance  wave- 
length. The  relative  growth  of  disturbance  amplitude  A,  defined  in 
Section  IID  Is  plotted  versus  time  for  each  case  on  Fig.  7.  We  see 
that  the  case  with  fewer  cells  grows  somewhat  slower  than  the  more 
accurate  case,  but  it  still  shows  the  same  basic  trend.  We  may  there- 
fore use  fewer  cells  when  we  are  Interested  in  the  trends  of  a series 
of  different  problems.  This  will  save  computer  time  and  enable  us  to 
study  the  results  much  sooner. 

5.  The  Effect  of  Disturbance  Amplitude 

In  our  previous  calculations,  for  the  most  part,  the  value  of 
the  amplitude  of  the  Initial  disturbance  has  been  taken  as  5%  of  the 
initial  jet  radius.  We  next  examine  the  case  where  the  Initial  distur- 
bance amplitude  Is  much  smaller,  O.IZ  of  the  Initial  radius.  Figure  8 
shows  the  results  of  a comparison  between  two  Identical  calculations 
except  for  the  amplitude  of  the  disturbance.  In  Fig.  8,  the  minimum 
radius  versus  time  Is  plotted  for  each  case.  We  see  that  after  a given 
amount  of  time,  the  case  with  the  smaller  Initial  amplitude  does  not 
neck  down  as  much  as  the  case  with  larger  Initial  amplitude,  but  signif- 
icant growth  does  occur.  In  fact,  percentage  wise,  the  growth  for  the 
smaller  amplitude  disturbance  Is  somewhat  larger  as  can  be  seen  by 
computing  the  relative  growth  A.  For  example,  at  t ■ 24  ps,  A for  the 
case  with  small  amplitude  Is  4.3,  whereas  A for  the  case  with  large 
amplitude  Is  3.6.  Recall  that  for  the  example  In  Section  IID  there  was 
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Comparison  of  HEMP  calculations  of  a stretching  jet  with  two  different 
mesh  sizes. 


t(ys) 

Figure  8.  A comparison  of  minimum  jet  radius  histories  as  predicted  by  HEMP  for 
t%K)  Jets  with  different  initial  disturbance  amplitudes. 


very  little  difference  in  A for  the  two  amplitudes  used.  Thus,  we  see 
that  even  for  this  small  amplitude  a definite  stability  phenomenon 
occurs.  We  also  should  note  that  for  the  times  considered,  all  the 
cells  remain  plastic  In  the  small  amplitude  case.  For  the  larger  ampli- 
tude case,  an  elastic  region  Is  observed  to  begin  at  approximately 
t ■ 12  us. 

In  our  previous  work  [4] , we  studied  the  effect  of  wavelength  on 
disturbance  growth  with  the  larger  amplitude.  In  the  light  of  the  results 
Just  presented,  we  conducted  a similar  wavelength  study  with  a smaller 
disturbance  amplitude.  A Jet  element  having  an  Initial  radius  of  1.5mm 
and  an  initial  strain  rate  of  0.029  us~^  was  used.  The  disturbance 
amplitude  was  0.1%  of  the  Initial  radius.  Four  separate  surface  distur- 
bances having  Initial  wavelengths  of  2.25mm,  4.5mm,  9mm,  and  18mm  were 
Introduced  and  HEMP  calculations  made  for  each  case.  For  convenience,  a 
smaller  number  of  cells  were  used  In  this  comparison  (2  cells  across  the 
radius).  Note  that  we  have  chosen  the  problem  such  that  the  case 
X « 4.5mm  will  eventually  grow  to  the  observed  segment  length  after  the 
experimentally  determined  breakup  time  for  a typical  shaped  charge  Jet. 

Let  us  denote  this  "correct"  Initial  wavelength  as  X “ 4.5mm.  We  can 
then  denote  all  the  cases  as  X • Xj./2,  X » X , X » 2a^,  and  X - 
The  relative  growth,  as  predicted  by  the  HEMP  code  for  each  of  these 
cases.  Is  plotted  vs.  time  In  Fig.  9.  From  this  plot  we  see  that  the 
case  X - Xj,  grows  significantly  faster  than  the  other  three  cases.  Thus, 
even  In  this  small  amplitude  case  where  the  material  remains  plastic 
throughout  and  no  elastic  unloading  occurs,  there  exists  a very  definite 
preferred  or  critical  disturbance  wavelength.  We  therefore  conclude 
that  the  critical  wavelength  Is  not  governed  by  elastic  unloading  In  the 
stretching  Jet. 


6.  The  Effect  of  Elasticity 

In  the  previous  section  we  have  Indicated  that  the  critical 
wavelength  of  disturbance  growth  In  shaped  charge  Jets  Is  a phenomenon 
which  does  not  depend  on  the  elastic  unloading  which  occurs  In  ttie  Jet. 
To  observe  the  effects  of  elasticity  further,  we  have  made  a study  of 
Jets  which  are  Identical  In  all  respects  except  for  the  elastic  moduli. 

In  HEMP,  elastic  moduli  which  are  required  as  Input  are  the  shear 
modulus  G and  the  bulk  modulus  K.  Note  that  K Is  entered  through  the 
equation  of  state  of  the  material.  The  various  combinations  of  elastic 
moduli  used  are  shown  In  Table  III. 
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Table  III 


Combinations  of  Elastic  Moduli  Used  In  Calculations 


Case  G 

No.  (102g  Pa) 

K 

(lO^G  Pa) 

E 

(lO^G  Pa) 

V 

1 0.456 

1.190 

1.213 

0.330 

2 0.456 

0.334 

0.940 

0.031 

3 0.456 

11.400 

1.350 

0.480 

4 0.410 

9.752 

1.213 

0.479 

5 0.590 

0.428 

1.213 

0.028 

6 0.225 

0.600 

0.600 

0.333 

7 0.750 

2.000 

2.000 

0.333 

Ncte:  1 M bar  - lO^G  Pa 

For  convenience,  Young's  modulus  E and  Poisson's  ratio  v are  also  listed 

Case  1 gives  the  standard  elastic  moduli 

of  copper 

we  have  been  using. 

In  Cases  2 and  3,  G Is  fixed  at  the  standard  value 

and  we  chose  two 

values  for  E noting  that  there  are  practical  bounds 

on  E given  by 

2G  < E < 3G.  These  bounds 

on  E come  directly  from 

the  practical  bounds 

on  Poisson's  ratio  0 < v < 

*j*.  In  Cases  4 and  5, 

we  have  fixed  E 

at  the  standard  value  and  let  G vary  between  the  practical  bounds 

(E/3)  < G < (E/2).  Finally 

, to  obtain  a 

larger  variation  In  G and  E we 

allowed  both  to  vary  keeping  v fixed. 

For  each  case,  the  Initial  configuration  shown 

In  Fig.  7 was 

used.  The  larger  cell  size 

(2  X 6)  was  also  used  In  all  cases.  The 

results  have  shown  that  the 

variations  In  elastic  moduli  do  not  serl- 

ously  affect  the  growth  of 

the  surface  disturbance. 

This  can  be  seen 

In  Table  IV,  where  we  have 

listed  the  relative  growth  and  the  minimum 

radius  at  t • 28  ys  for  each  of  the  cases  studied. 

Note  that  the 

largest  deviation  In  A from 

the  standard 

value  Is  5%. 

Table  IV 

Results 

of  Elastic  Moduli  Study 

t ■ 28  ys 

Case 

'^mln 

A 

No. 

(mm) 

1 

.672 

3.81 

2 

.663 

3.89 

3 

.652 

4.01 

4 

.647 

4.06 

5 

.663 

3.88 

6 

.663 

3.90 

7 

.676 

3.77 

* Note  theoretically  -H  < 

V < 4;  however,  no  real 

materials  have  been 

found  with  V < 0,  therefore  we  Impose 

practical 

bounds  of  0 < V < 

30 
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III.  TWO-DIMENSIONAL  CALCULATIONS  OF  MULTIPLE  INTERIOR  IMPACT 

PENETRATION 

Another  topic  examined  under  the  present  study  Is  the  penetration 
of  broken  shaped  charge  jets.  The  object  of  this  portion  of  the  study 
Is  to  define  the  mechanisms  which  cause  the  decrease  In  penetration 
associated  with  the  Increase  In  the  gap  between  jet  particles  as  exper- 
imentally observed  In  [8]. 

A.  Background 

Under  a recent  program  conducted  at  Drexel  University  [7-10], 
experiments  were  performed  which  studied  the  effects  of  projectiles 
Impacting  upon  the  bottom  of  a pre-drllled  hole  in  a target.  We  define 
the  term  "multiple  Interior  Impact"  to  describe  this  general  setup. 

These  tests  were  designed  to  simulate  penetration  by  a segmented  shaped 
charge  jet.  In  these  experiments,  two  low  speed  (approximately  Imm/ps) 
lead  projectiles,  separated  by  a pre-determined  gap,  were  Impacted  upon 
soft  homogeneous  lead  targets.  It  was  concluded  from  these  experiments 
that  a definite  relationship  exists  between  the  size  of  the  gap  and 
total  penetration.  As  the  gap  between  projectiles  increases,  total 
penetration,  in  general,  decreases.  The  probable  cause  for  this  decrease 
was  attributed  to  the  Interference  between  the  trailing  projectile  and 
the  occlusion  or  "back  flow"  caused  by  the  leading  projectile.  This  In- 
terference causes  the  formation  of  three  distinct  cavities  in  the  target 
as  shown  In  [8-10].  It  is  likely  that  this  Interference  mechanism  is 
the  cause  of  the  "scalloped"  hole  observed  from  the  penetration  of 
broken  shaped  charge  jets  (see  [29]). 

During  the  present  study  we  conducted  two-dimensional  HELP  code 
calculations  to  gain  an  understanding  of  the  multiple  Interior  impact 
process.  In  our  first  series  of  calculations  we  used  lead  projectiles 
at  low  speeds,  impacting  lead  targets  to  simulate  experiments  in  [8-10]. 

We  also  computed  copper  projectiles,  having  dimensions  of  typical 
shaped  charge  jet  segments.  Impacting  at  typical  shaped  charge  jet 
velocities  onto  steel  targets. 

B.  General  Approach 

The  HELP  code  [11]  was  used  for  all  the  multiple  Interior  Impact 
calculations.  This  two-dimensional  finite-difference  code  uses  a com- 
bined Eulerlan-Lagrangian  hydrodynamic  method  which  treats  multi- 
material  problems  with  a continuous  description  of  the  moving  mass. 

The  method  used  In  HELP  handles  material  Interfaces  and  free  surfaces 
by  propagating  these  surfaces  as  sharp  Lagranglan  discontinuities 
through  the  stationary  Eulerlan  grid.  This  code  has  the  capability  of 
treating  elastic  and  plastic  effects. 

The  version  of  HELP  used  did  not  have  the  necessary  routines  to  set 
up  the  target  and  projectile  shapes  used  In  the  Drexel  experiments  [8-10]. 
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Therefore  an  auxiliary  set  up  package,  knovm  as  HELPCLAM  (30],  obtained 
through  the  courtesy  of  John  Harrison,  BRL,  was  used  to  set  up  the 
necessary  geometries.  The  original  projectiles  were  modeled  as  cylin- 
ders or  cylinders  with  conical  tips.  An  example  of  the  computational 
grid  and  a typical  Initial  configuration  used  for  these  problems  are 
shown  in  Figures  10  and  11  respectively.  Table  V summarizes  the 
Initial  conditions  for  the  four  configurations  studied. 

The  Tlllotson  equation  of  state,  as  specified  In  the  HELP  code 
(10)  for  each  of  the  materials  considered,  was  used. 

C.  Impact  of  Lead  Projectiles  on  Lead  Targets 

We  have  computed  three  separate  configurations  of  low  speed 
( * Imm/psec)  multiple  interior  Impact  problems  using  projectiles  and 
targets  made  of  lead.  The  first  calculation  (Table  V,  No.  1),  simu- 
lated one  of  the  experiments  presented  In  [9].  Figure  11  shows  the 
Initial  configuration  for  this  problem.  The  Initial  gap  between  pro- 
jectiles Is  approximately  22mm.  The  yield  strength  for  lead  was  chosen 
to  be  6.89  MPa. 


Table  V 

Summary  of  the  Initial  conditions  of  the  HELP  Code  Calculations 


Run  No. 

Projectile 

Projectile 
yield  strength 
(Pa) 

Gap 

(mm) 

Target 

Target 

yield  strength 
(Pa) 

1 

Lead 

6.89x10^ 

22 

Lead 

6.89x10^ 

2 

Lead 

6.89x10^ 

37 

Lead 

2.756x10^ 

3 

Lead 

6.89x10^ 

37 

Lead 

2.756x105 

4 

Copper 

6.89x10^ 

14 

Steel 

2.41x10® 

Figures  12  through  16  show  the  computed  projectile  penetration 
at  successive  times  after  Impact.  Notice  that  these  computer  simula- 
tions show  no  Interference  between  the  trailing  projectile  and  the 
target  "back  flow"  caused  by  the  first  projectile.  Thus,  the  computer 
simulation  does  not  show  the  formation  of  the  three  distinct  cavities 
observed  in  the  experiments  of  [9], 

A plot  of  computed  depth  of  penetration  vs.  time  after  impact  is 
given  in  Figure  17.  The  total  penetration  depth  measured  for  the 
comparable  experiment  [9]  is  approximately  42mm.  As  one  can  see  the 
HELP  predicted  penetration  is  only  about  one-half  of  the  experimental 
value. 

Though  no  interference  was  predicted  in  these  calculations,  the 
portion  of  the  target  directly  parallel  to  the  target  hole  axis  moves 
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Itiple  interior  impact  problems. 


Axis  of  Symmetry 


Figure  11.  Initial  configuration  of  the  interior  impact  of  lead  project 
onto  a lead  target. 


Axis  of  Symmetry 


Figure  13.  Configuration  of  lead  interior  impact  at  30ps  after  initial  impact  (initial  gap 


Configuration  of  lead  interior  Impact  at  A7gs  after  initial 


Axis  of  Symmetry 


Figure  16.  Configuration  of  lead  interior  Impact  at  71ps  after  initial  Impact  (Initial  gap  = 22mm). 


radially  inward  as  can  be  seen  in  Figures  12  through  16.  If  the  gap 
between  projectiles  were  Increased,  this  region  of  the  target  would 
have  more  time  to  move  radially  Inward  to  a position  where  it  would 
Interfere  with  the  path  of  the  trailing  projectile.  Therefore,  a 
second  run  was  made,  this  time  modeling  an  experiment  with  a larger 
spacing  between  projectiles. 

This  calculation  (Table  V,  No. 2)  simulates  one  of  the  experiments 
reported  in  [8].  The  initial  gap  between  projectiles  was  approxi- 
mately 37mm.  From  the  discrepancy  observed  in  the  total  penetration 
of  the  previous  case,  we  felt  that  the  lead  was  not  being  modelled 
properly  in  the  code.  The  yield  strength  of  the  lead  target  was  thus 
lowered  to  2.76  x 10^  Pa  in  an  attempt  to  alleviate  this  problem. 

These  and  all  other  Important  dimensions  are  given  in  Figure  18, 

Figures  19  through  21  show  the  calculated  penetration  at  succes- 
sive times  after  impact.  From  these  computer  simulations  we  have 
calculated  that  there  is  not  enough  time  for  the  target  "back  flow"  to 
move  radially  inward  and  Interfere  with  the  trailing  projectile. 

In  the  next  run  (Table  V,  No. 3),  the  target  yield  strength  was 
lowered  by  a factor  of  ten,  to  2.76  x 10^  Pa  and  computations  were 
made.  All  other  initial  conditions  were  the  same,  as  indicated  in 
Figure  18,  Figures  22  through  24  are  the  computer  simulations  of  this 
run  at  successive  times  after  Impact.  Again  the  radial  velocity  cal- 
culations Indicate  that  no  Interference  will  occur  between  the  trail- 
ing projectile  and  the  target. 

At  24  psec  after  the  initial  Impact  the  calculated  penetration 
was  10.38  mm  as  compared  to  9.42  mm  penetration  for  the  previous  case 
(Table  V,  No. 2).  This  Increase  in  penetration  is  not  substantial  con- 
sidering the  90Z  decrease  in  target  yield  strength. 

The  HELP  code  calculations  did  not  compare  well  with  the  experi- 
mental results  of  [8]  and  [9].  The  computed  penetration  was  only  about 
one-half  of  that  measured.  The  code  did  not  predict  interference  be- 
tween the  trailing  projectile  and  the  occlusion  caused  by  the  impact  of 
the  first,  and  thus  the  subsequent  triple  cavity  in  the  target  was  not 
observed.  One  possible  explanation  for  these  computations  being  un- 
successful is  the  anomalous  behavior  of  lead  or,  more  appropriately, 
the  lead  equation  of  state  used  in  the  code. 

D.  Impact  of  Copper  Prolectlles  on  Steel  Targets 

In  an  attempt  to  understand  the  mechanism  of  this  gap  effect  be- 
tween projectiles,  the  HELP  code  was  used  to  calculate  the  interior 
impact  of  two  copper  projectiles  on  a steel  target.  The  dimensions  of 
the  copper  projectiles  were  chosen  to  simulate  typical  jet  segments  of 
an  81.3  mm  BRL  charge.  These  dimensions,  along  with  those  of  the  tar- 
get, are  given  in  Figure  25.  Figures  26  and  27  are  the  computer  simu- 
lations of  this  run  at  4 and  9 psec  after  impact  respectively. 
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Axis  of  Symmetry 


lead  target. 


Axis  of  Symmetry 


Figure  19.  Configuration  of  the  lead  Interior  Impact  (Initial  gap  = 37mm)  at  8us 
after  Initial  Impact  as  predicted  by  HELP  code. 
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Figure  20.  Configuration  of  the  lead  interior  Impact  (initial  gap  = llmm)  at  13ys 
after  initial  Impact. 


Axis  of  Symmetry 


Figure  22.  Configuration  of  the  lead  interior  impact  (Initial  gap  • 37nim)  with  a lower 

target  yield  strength,  at  19vs  after  initial  impact  as  predicted  by  HELP  code. 


Figure  23.  Configuration  of  lead  interior  impact  (initial  gap  = 37iran,  lower  yield  strength) 
at  23us  after  initial  impact. 
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at  29iis  after  initial  Impact. 


SCALE 


Figure  25.  Initial  configuration  of  the  interior  impact  of  copper  projectiles  onto  a steel  target 
(initial  gap  = 14mm). 


In  thli  set  of  computations  we  successfully  showed  that  the  trail- 
ing projectile  does  Indeed  Interfere  with  the  occlusion  or  "back  flow" 
caused  by  the  Impact  of  the  first.  In  Figure  26  the  trailing  projectile 
has  already  Interfered  with  the  occlusion  and  that  portion  of  the  target 
has  begun  to  move  radially  outward.  The  actual  Interference  occurred  at 
a time  prior  to  4 ps.  Note  that  In  an  Eulerlan  code,  once  two  materials 
enter  the  same  cell,  the  Interfaces  within  that  cell  are  assumed  to 
impact  or  be  In  contact.  This  explains  why.  In  Figure  26,  at  first 
glance,  one  might  conclude  that  no  Impact  occurred  between  the  trailing 
projectile  and  target  occlusion  but,  upon  further  Inspection  of  the  de- 
tailed calculations,  the  radial  movement  of  the  target  outward,  due  to 
this  Impact,  becomes  clear.  The  formation  of  the  three  distinct  cavities 
In  the  target  can  be  observed  In  Figure  27.  This  "scalloped"  effect  Is 
observed  In  penetration  tests  of  broken  shaped  charge  jets  [29]. 
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IV.  CONCLUSIONS 


We  will  now  briefly  present  the  major  conclusions  on  the  two  tasks 
of  this  study,  jet  breakup  and  broken  jet  penetration.  We  will  also 
give  recommendations  for  continued  work  on  these  topics. 

A.  Stability  and  Breakup  of  Jets 

In  the  study  of  jet  breakup,  we  have  formulated  an  analytical 
approach  and  have  continued  the  numerical  computations  of  our  previous 
study  [4].  A one-dlmenslonal  theory  was  formulated  to  solve  the  stretch- 
ing jet  problem  subject  to  surface  disturbances.  The  theory  uses  a 
strain  hardening  plastic  model  and  assumes  the  surface  disturbances  to 
be  small.  Elastic  loading  and  unloading  is  not  Included  in  the  develop- 
ment. The  theory  yields  a stability  criterion  which  depends  on  the 
stress-strain  curve.  For  an  unstable  jet,  any  small  disturbance  will 
grow.  The  theory  also  predicts  the  existence  of  a critical  distur- 
bance wavelength.  That  Is,  disturbances  of  this  wavelength  will  grow 
faster  than  any  other  for  a given  set  of  conditions.  For  disturbances 
having  this  critical  wavelength,  the  stability  criterion  states  that  at 
those  times  when  the  true  stress  Is  smaller  than  or  equal  to  the  slope 
of  the  true  stress-true  strain  curve,  the  jet  Is  stable;  otherwise.  It 
is  unstable.  For  a perfectly  plastic  jet,  the  critical  wavelength  Is 
always  given  by  ^^g^^’^O^crlt  “ ^ ^ regardless  of  material  properties 

or  strain  rate.  Furthermore,  for  perfectly  plastic  jets.  If  A^/Zr^  > ti/Z, 

the  jet  Is  always  unstable.  The  theory  also  immediately  reveals  that 
disturbance  growth  rate  in  the  jet  Increases  as  the  strain  rate  no 
decreases.  Finally,  the  disturbance  growth  Increases  with  the  ratio  of 
flow  stress  to  density,  Oq/Pq. 

Comparison  with  two-dimensional  calculations  Indicates  that  the 
linearized  one-dlmenslonal  theory  predicts  all  the  pertinent  trends  In- 
dicated by  the  two-dimensional  results.  The  jet  dimensions  and  strain 
rates  used  In  the  calculations  are  comparable  to  those  observed  In  a 
typical  jet  from  a BRL  81.3mm,  42®,  copper-lined  conical  charge.  The 
growth  calculated  from  the  theory  appears  reasonably  accurate  when  the 
disturbance  Is  Initiated  In  the  later  stages  of  the  shaped  charge  jet 
elongation  process  (l.e.  when  the  strain  rate  Is  relatively  low).  For 
the  earlier  times  after  jet  formation  (l.e.  when  the  strain  rate  Is 
relatively  high)  the  theoretical  results  are  not  adequate.  Finally 
we  note  that  the  two-dimensional  results  Indicate  a dependence  of 
critical  wavelength  to  diameter  ratio  on  the  Initial  strain  rate.  We 
have  found  Xq/2tq  to  vary  between  1.5  for  higher  strain  rates  to  2.2 
for  lower  strain  rates  In  the  jets  studied  here.  However,  the  linear- 
ized theory  predicts  a value  of  Xq/2tq  - 2.2  regardless  of  strain  rate. 

At  present,  we  feel  that  these  problems  can  be  resolved  by  providing  a 
more  accurate  solution  to  the  theory.  One  method  which  appears  very 
promising  Is  to  apply  the  linearized  solution  In  a stepwise  fashion  for 
small  Increments  of  time  taking  care  to  match  the  proper  conditions  at 
each  step. 
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Tntt  fcll->wlng  r.ce  recoaoended  areaa  for  future  work  on  Jet  atablllty 
and  breakup: 

1.  The  theory  ahould  be  Improved  to  account  for  growth  at 
longer  times  after  disturbance  Initiation,  higher  strain  rates, 
and  the  affect  of  strain  rate  on  critical  wavelengtha,  as 
mentioned  above,  one  approach  would  be  to  apply  the  present 
theory  In  a stepwise  fashion  for  small  Increments  of  time  and 
match  the  proper  conditions  between  steps. 

2.  If  the  Initial  disturbance  amplitude  and  Its  time  of 
Initiation  are  known,  we  can  calculate  the  wavelength  and 
amplitude  at  any  later  time.  However,  In  a real  Jetting  sltu" 
atlon,  we  do  not  know  precisely  when  to  Initiate  the  dlst'ir- 
bance  and  at  what  amplitude.  This  must  be  know^  *'o  predict 
the  proper  number  of  eventual  Jet  particles. 

3.  We  currently  do  not  have  a reliable  method  to  determine 
how  large  the  amplitude  of  a disturbance  must  grow  to  cause 

a break.  "Brittle"  Jets  and  "ductile"  Jets  may  behave  differently. 
This  Information  Is  necessary  to  predict  the  eventual  aspect 
ratio  of  the  Jet  segments  and  the  precise  breakup  time. 

4.  At  present,  we  do  not  know  precisely  the  Instantaneous  Jet 
properties  under  Impulsive  and  dynamic  stretching  conditions. 
Yielding,  work  hardening,  rate  effects,  etc.  will  all  play  an 
Important  part  In  predicting  disturbance  growth  and  breakup  in 
actual  shaped  charge  Jets. 

B.  Penetration  Calculations 


Various  calculations  of  multiple  Interior  impact  were  conducted 
to  determine  the  isechanlsms  Involved  In  the  penetration  of  broken  Jets. 
The  calculations  Involving  lead  projectiles  and  targets  did  not  show 
the  Interference  phenomenon  expected.  In  fact,  calculated  penetrations 
were  approximately  50^  smaller  than  chose  observed  in  the  experiments 
of  [9].  At  present.  It  Is  felt  that  this  discrepancy  Is  due  to  the 
anomalous  behavior  of  the  lead  equation  of  state  used  In  the  code. 

In  the  calculations  of  copper  projectiles  onto  steel  targets,  Che 
expected  Interference  phenomenon  was  observed  and  the  formation  of 
three  distinct  cavities  from  two  projectiles  was  shown.  This  may  ex- 
plain the  "scalloped"  hole  formed  when  a target  Is  penetrated  by  a 
broken  shaped  charge  Jet. 

Though  Che  expected  penetration  phenomenon  was  observed  In  the 
copper-steel  case,  the  calculated  hole  size  was  smaller  than  expected. 
An  attempt  will  be  made  to  find  out  why  this  discrepancy  exists  and. 

If  successful,  any  necessary  modifications  of  Che  code  will  be  made  so 
Chat  In  future  calculations  It  will  be  able  to  more  accurately  model 
these  Impact  problems. 


The  next  step  is  to  compute  the  same  copper-steel  problem  with  no 
gap  between  the  projectiles.  This  will  demonstrate  whether  the  pene- 
tration does  Indeed  decrease  significantly  with  increase  In  gap  as  shown 
In  [7]  and  [8],  If  it  turns  out  that  this  gap  effect  Is  significant, 
then  current  penetration  theories,  which  add  up  the  segment  length  to 
compute  penetration,  should  be  modified  to  Include  a factor  which  also 
accounts  for  the  gap. 

Finally,  It  Is  recommended  that  future  calculations  Include  lower 
speed  Impacts  (2-A  km/sec)  of  copper  projectiles  on  steel  targets  aimed 
at  determining  the  proper  "cutoff"  criterion  of  jet  penetration.  Also 
HELP  code  penetration  calculations  of  continuous  jets  having  a velocity 
gradient  along  their  length  should  also  be  undertaken. 
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Appendix 


In  this  appendix  we  will  present  the  stress  distribution  for 
typical  stretching  Jets  subjected  to  a surface  disturbance.  Two 
stretching  Jet  segments  having  a slight  difference  in  initial  dis- 
turbance and  a slight  difference  in  the  initial  mesh  were  simulated 
using  the  HEMP  code.  The  stresses  calculated  are  presented  in  both 
tabular  and  graphical  form. 
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A.  Tabular  Display  of  Stresses  as  Predicted  by  HEMP  Code  for  a 
Typical  Stretching  Jet. 


In  this  section  we  tabulate  stress  data  as  computed  by  HEMP  for  the 
problem  depicted  In  Fig.  Al.  The  Initial  data  for  this  problem  Is  as 
follows; 


3 

Initial  density  Pq  “ 8.9  gm/cm  yield  strength  - 0.2  GPa 
Initial  pressure  P^  - 0 shear  modulus  G - 45.6  GPa 


Figure  Al . Initial  geometry  and  grid  layout  for  a typical  HEMP  calculation. 


For  this  problem,  the  stresses  as  predicted  by  HEMP  are  listed  as 
follows: 
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B.  Graphical  Presentation  of  Axial  Stresses  In  a Typical 
Stretching  Jet 


In  this  section  we  will  graphically  display  the  stress  data  as 
calculated  using  the  HEMP  code  for  the  problem  depicted  In  Figure  A2. 
All  other  Initial  data  the  same  as  the  previous  problem  In  part  A of 
appendix. 

The  following  pages  contain  plots  of  axial  stress  (a^)  vs.  time 
for  the  cells  specified. 


Figure  A2.  Initial  geometry  and  grid  layout  for  a typical  HEMP  calculation 
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For  cell  (6,5) 


Figure  A3.  Axial  stress  vs.  time  plots  for  cells  (6,5)  and  (8,5)  as  defined  in  Figure  A2 


in  Tigure 


For  cell  (6,3) 


M 


Figure  A5.  Axial  stress  vs.  time  plots  for  cells  (6,3)  and  (8,3)  as  defined  in  Figure  A2 


(1,3) 


Figure  A6.  Axial  stress  vs.  time  plots  for  lells  (1,3)  and  (4,3)  as  defined  in  Figure  A2 . 


DISTRIBUTION  LIST 


No,  ot  No.  of 

Copies  Orj^anization  Copies  Organization 


12  Commander  2 

Defense  Documentation  Center 
ATTN : DDC-TCA 
Cameron  Station 
Alexandria,  VA  22314 

I Commander 

US  Array  Materiel  Development  I 
and  Readiness  Command 
ATTN:  DRCDMA-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

6 

1 Commander 

US  Army  Aviation  Research 
and  Development  Comn and 
ATTN ; DRSAV-E 
12th  and  Spruce  Streets 
St.  Louis,  MO  63166 

1 Director 

US  Army  Air  Mobility  Research 
and  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035  1 

1 Commander 

US  Army  Electronics  Command 

ATTN:  DRSEL-RD 

Fort  Monmouth,  NJ  07703 

2 

1 Commander 

US  Army  Missile  Research 
and  Development  Command 
ATTN:  DRDMI-R 
Redstone  Arsenal,  AL  35809 

1 Commander  1 

US  Army  Tank  Automotive 
Development  Command 
ATTN:  DRDTA-RWL 
Warren,  MI  48090 


Commander 

US  Army  Mobility  Equipment 

Research  8 Development  Command 
ATTN:  Tech  Docu  Cen,  Bldg.  315 
DRSML-RZT 

Fort  Belvoir,  VA  22060 
Commander 

US  Army  Armament  Materiel 
Readiness  Command 
Rock  Island,  IL  61202 

Commander 

US  Army  Armament  Research 
and  Development  Command 
AITN:  DRDAR 

Mr.  T.  Stevens 
Mr.  G. Randers-Fehrson 
Dr.  N.  Clark 
Mr.  J.  Mershkowitz 
Mr.  J.  Pearson 
Tech  Lib 
Dover,  NJ  07801 

Commander 

US  Army  Harry  Diamond  Labs 
ATTN : DRXDO-TI 
2800  Powder  Mill  Road 
Adelphi,  MD  20783 

Commander 

US  Army  Materials  and 

Mechanics  Research  Center 
ATTN:  DRXMR-RF,  J.  Mescall 
Tech  Lib 

Watertown,  MA  02172 
Director 

US  Army  TRADOC  Systems 
Analysis  Activity 
ATTN:  ATAA-SL  (lech  Lib) 

White  Sands  Missile  Range 
NM  88002 


69 


DISTRIBUTION  LIST 


No.  of 

Copies  Or^ani zat ion 

1 Assistant  Secretary  of 

the  Army  (R^D) 

ATTN:  Asst  for  Research 
Washington,  DC  20310 

2 HQDA  (UAMA-ZA;  DANIA-AR) 
Washington,  DC  20310 

1 Commander 
US  Army  Research  Office 
P.  0.  Box  12211 
Research  Triangle  Park 
NC  27709 

2 Commander 
US  Naval  Air  Systems  Command 
ATTN:  Code  AIR- 3 10 

Code  AIR-350 
Washington,  DC  203b0 

1 Commander 

US  Naval  Ordnance  Sys  Command 
ATTN:  Code  ORU-0332 
Washington,  DC  20360 

2 Chief  of  Naval  Research 
Department  of  the  Navy 
ATTN:  Code  427 

Code  470 

Washington,  DC  20325 


No.  of 

Copies  Organization 

1 Commander 

US  Naval  Weapons  Center 
ATTN:  Code  45,  Tech  Lib 
China  Lake,  CA  93555 

1 Commander 

US  Naval  Research  Laboratory 
Washington,  DC  20375 

I USAF/APRDDA 

Washington,  DC  20311 

1 AFSC/SDW 
Andrews  AFB 
Washington,  DC  20311 

1 US  Air  Force  Academy 
ATTN:  Code  F.JS-41  (NC) 

Tech  Lib 

Colorado  Springs,  CO  80840 

I AFWL  (SUL,  LT  Tennant) 

Kirtland  AFB,  NM  87116 

1 AFLC/MMWMC 

Wright-Patterson  AFB,  011  45433 
1 AFAL/AVIV 

Wright-Patterson  AFB,  OH  45433 
6 Director 

Lawrence  Livermore  Laboratory 
ATTN:  Dr.  .1.  Kury 

Dr.  M.  Wilkins 
Dr.  E.  Lee 
Dr.  U.  Ilornig 
Dr.  .).  Knowles 
Tech  Lib 
P.  0.  Box  808 
Los  Alamos,  NM  87544 


1 Commander 

US  Naval  Surface  Weapons  Center 
ATTN:  Code  730,  Lib 
Silver  Spring,  MD  20910 

2 Commander 

US  Naval  Surface  Weapons  Center 
ATTN:  Code  DG-50 
Tech  Lib 

Dahlgren,  VA  22448 


70 


0 


t 


DISTRIBUTION  LIST 


No.  of 
Copies 


Organization 


No.  of 
Copies 


Organization 


4 


1 


Battelle 

Columbus  Laboratories 
ATTN:  Mr.  Joseph  L.  Backofen 
505  King  Avenue 
Columbus,  on  43201 

Dyna  East  Corporation 
ATTN:  P.  C.  Chou 
J.  Carleone 
C.  A.  Tanzio 
R.  D.  Cicarelli 
227  Hemlock  Road 
Wynnewood,  PA  1909b 

Honeywell,  Inc. 

Government  and  Aeronautical 
Products  Division 
ATTN:  C.  R.  Hargreaves 
600  Second  Street 
Hopkins,  MN  55343 

Physics  International  Corp 
ATTN:  Dr.  L.  Berhman 
2700  Merced  Street 
San  Leandro,  CA  94577 


1 Systems,  Science  5 Software 
ATTN:  Ur.  R.  Sedgwick 
P.  0.  Box  1620 
La  Jolla,  CA  92037 

1 Drexel  Institute  of  Technology 
Wave  Propagation  Research 

Center 

ATTN:  Prof.  P.  Chou 
Philadelphia,  PA  19104 

2 University  of  California 
Los  Alamos  Scientific  Lab 
ATTN:  Dr.  J.  Walsh 

Tech  Lib 
P.  0.  Box  1663 
Los  Alamos,  NM  87544 

1 University  of  Dayton  Research 
Institute 

ATTN:  Ur.  Stephen  J.  Bless 
Dayton,  OH  45469 


1 Sandia  Laboratories  Marine  Corps  Ln  Ofc 

ATTN:  Dr.  W.  Herrmann  Dir,  US/\MSAA 

Albuquerque,  NM  87115 


1 Shock  Hydrodynamics 
ATTN:  Dr.  L.  Zernow 
4710-4716  Vineland  Avenue 
North  Hollywood,  CA  91602 


